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We use analytic calculations of the post-recombination gravitational effects of cosmic strings to 
estimate the resulting CMB power spectrum, bispectrum and trispectrum. We place a particu- 
lar emphasis on multipole regimes relevant for forthcoming CMB experiments, notably the Planck 
satellite. These calculations use a fiat sky approximation, generalising previous work by integrat- 
ing string contributions from last scattering to the present day, finding the dominant contributions 
to the correlators for multipoles / > 50. We find a well-behaved shape for the string bispectrum 
(without divergences) which is easily distinguishable from inflationary bispectra which possess sig- 
nificant acoustic peaks. We estimate that the nonlinearity parameter characterising the bispectrum 
is approximately }nl ~ —20 (given present string constraints from the CMB power spectrum. We 
also apply these unequal time correlator methods to calculate the trispectrum for parrallelogram 
configurations, again valid over a large range of angular scales relevant for WMAP and Planck, as 
well as on very small angular scales. We find that, unlike the bispectrum which is suppressed by 
symmetry considerations, the trispectrum for cosmic strings is large. Our current estimate for the 
trispectrum parameter is tnl ~ 10 4 , which may provide one of the strongest constraints on the 
string model as observational estimates for the trispectrum improve. 

INTRODUCTION 

Cosmic strings are a common feature in fundamental cosmology scenarios, such as brane inflation [1J, and they 
appear to be generic in realistic grand unified theories [2J. Cosmic strings leave a distinct 'line-like' signature in the 
cosmic microwave background (CMB) [3J H] which offers one of the best prospects for their detection. Strings create 
this imprint after recombination by perturbing photons through relativistic gravitational effects. While inflationary 
fluctuations are believed to dominate the overall CMB power spectrum, current constraints suggest that up to 10% 
of the signal could be contributed by cosmic strings [SI E] • A higher proportion is incompatible with WMAP because 
the cosmic string power spectrum around multipoles I ~ 200 is dominated by metric fluctuations with a relatively 
featureless spectrum, rather than the pre-recombination acoustic peaks characterising inflation. The gravitational 
strength of cosmic strings is determined by the parameter Gli = (r//mpi) 2 which is the ratio of the string tension \x 
to the square of the Planck mass mpi. The present WMAP limit on the string contribution to the CMB translates 
into a strong constraint on the parameter G/j, < 2.5 x 10 -7 [5]. 

Inflationary CMB fluctuations in the standard picture are very nearly Gaussian, so higher order correlators offer the 
prospect of further differentiation between these and competing signals from cosmic strings. The spherical harmonic 
transform of the three-point CMB correlator is the bispectrum Bi ± i 2 i 3 . Usually discussions of the bispectrum are 
simplified further by focusing on the nonlinearity parameter fpjL which is roughly the ratio of the three point correlator 
to the square of the two-point correlator, that is, Jnl ~ (CCC)/(CC) 2 where the ( are the primordial curvature 
fluctuations (dimensionless) which seed the CMB anisotropics. For a perturbative theory like inflation, we expect 
the second-order fluctuations C*- 2 "* to be constructed from convolutions of the linear perturbations ~ C^ 1 * 1 * C*" 1 "*- 
This implies that the leading order term in the three-point correlator can be expected to behave as (C ^C^C^) ~ 
(C^C*' 1 ') 2 ' that is, /jvl ~ 1. In fact, for standard single field inflation there is considerable further suppression 
from slow-roll and the primordial signal /nl ~ O(0. 01) [7]. One can only obtain a significantly larger /jvl in more 
exotic inflationary models with multiple scalar fields or non-canonical kinetic terms (see, for example, the reviews in 
[5J[S]). Similar arguments apply to the trispectrum or four-point correlator characterised by tml ~ (CCCC)ACC) 3 with 
t nl 1 f° r standard single field inflation. 

In contrast, the cosmic string signature has an inherently non-perturbative origin so we do not expect Jnl to be 
small nor, more particularly, t^l- For the higher order correlators from cosmic strings (measured relative to the 
string power spectrum (£0cs)i the relevant scaling with the parameter G/j, is f^L ~ A(Gfj,)~ 1 |10j and, as we shall 
discuss here, tml ~ ^(Gli) 2 with coefficients A and B determined by geometric and dynamical considerations. There 
is a suppressed amplitude A for the bispectrum (due to symmetry cancellations [10 ]), and we contrast this here with 
one of the key results of this paper which is the much larger relative value of B for the trispectrum. Here, we shall 
compare the higher order correlators from strings we calculate with the dominant inflationary power spectrum (CC)m/ 
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and we will show that the trispectrum with tnl 3§> 1 is the more likely correlator to produce the strongest constraints 
on cosmic strings (or with which to identify them). 

In this paper, we use analytic calculations of the post-recombination string signature to estimate the power spec- 
trum, bispectrum and trispectrum relevant for forthcoming CMB experiments, particularly the Planck satellite. We 
generalise previous analytic calculations of Hindmarsh and collaborators for the power spectrum and bispectrum 
[TU] based on a flat sky approximation. While this was an important first step, the summation employed was essen- 
tially confined only to strings close to the surface of last scattering and it is only relevant for very small angular scales. 
Here, we note that the dominant contribution on large angular scales (as well as very small angles) comes from the 
same gravitational GKS effects but from strings at late times well after last scattering. By integrating in time over the 
unequal time correlators, we obtain an approximate bispectrum valid for all I > 50 (where the flat sky approximation 
breaks down), i.e. useful for both WMAP and Planck. We also apply these unequal time correlator methods to 
the first calculation of the trispectrum (for parralelogram configurations), which is again valid over a large range of 
angular scales. Our bispectrum calculation also differs from ref. |10j by eliminating the divergences they find in their 
Gaussian integrals for flattened triangles. We note that these are actually cut-off by the behaviour of other terms 
in the integrand (due to momentum constraints), yielding a finite result in this regime. The cross-sectional shape of 
the bispectrum we present is relatively featureless, except for a finite rise near the edges (for flattened triangles) and 
suppression towards the corners because of causality contraints (squeezed triangle limit) . This well-behaved shape is 
suitable for the bispectrum estimation methods being developed for Planck and is easily distinguishable from bispectra 
predicted by inflation because of the absence of acoustic peaks [8]. Of course, these analytic calculations neglect im- 
portant recombination effects which also provide significant contributions to the string bispectrum for 500 < I < 2000. 
However, determining the extent to which these contributions can confuse the non-Gaussian 'line-like' signature of 
cosmic strings will be the subject of future study [12"] . 

These results for the string bispectrum and trispectrum are important for CMB experiments and should be valid 
and dominant at both large and very small angular scales. While the Planck satellite does not have the resolution to 
see individual string signatures, it should be possible to obtain statistically significant constraints on cosmic strings 
that compete with limits on G/i from the power spectrum, especially for the trispectrum when these techniques are 
fully developed. On very small angular scales the string power spectrum begins to dominate over the inflationary 
signature (I > 3000) because it is not influenced by exponentially decaying transfer functions. Here direct detection 
of 'line-like' signatures may prove possible provided that experiments can achieve fiK sensitivities, as anticipated by 
AMI and ACT, for example. Again, searching in larger data sets for higher order correlators may provide statistically 
more significant results. 

GOTT-KAISER-STEBBINS EFFECT 

In this section we calculate the Gott-Kaiser-Stebbins (GKS) effect which is expected to give the principle contribu- 
tion to cosmic strings on subhorizon scales where we may use the flat sky approximation and can ignore the response 
of cosmological fluids. With these approximations the integrated Sachs Wolfe effect is equal to the temperature dis- 
continuity in the CMB. In fact we shall show in this paper that it is the 'late time small angular effect' induced by 
the GKS effect that also dominates large angular scales for cosmic strings. The discussion here otherwise follows that 

of[TU]. 

Consider a photon with 4- momentum = (E 1 , 0,0, E) and a string of coordinates X^(cr,t), where (a, t) are 
the worldsheet coordinates of the string (with the gauge chosen such that t corresponds to the (conformal) time 
coordinate). The unperturbed geodesies can be written as 

= + Xp" 

The perturbation to the energy along the photon path (which gives the ISW effect) is 

5 P ^ = - 1 -J\ up ^{z{\))pypd\ 

where h^ v is the metric perturbation which in general must be calculated. However, in refs 
out that we can simplify the discussion by considering 
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where V_l represents partial derivatives with respect to the transverse coordinates to the unperturbed geodesies of 
the photons. 

Now suppose that the photon motion is in the z-direction. This implies that (d t — d z )f(Z) = £? -1 ^. Also we note 
the following 

w 2 s Pp = (d 2 - d 2 - d 2 )s Pp 



= -I [ 1 dx[(d t + d z )(d t ~d z )-d 2 ]K p . p (z(x)) P » P ? 



Now suppose we use the harmonic gauge for the metric perturbations. Then we have d 2 h pv = 167rG(T piy — l/2ri fljl/ T) 
which implies that we have 



V 2 ^ 



+ 8nGd p f 1 d\T vp p v pP . 

i J An 
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Next we define p 11 = p^/E. Since we are interested in the energy we focus on the [i — component. We then obtain 1 

Iplfip — g ^ 1 PI,]PP 



T no" - — dT P l <°P l _ t &iP + T ■ ■ 



1 dTpi.pp 1 
9 iT E dX 



where we have used T p jj = T p0 fi in the second line and we define = d l —p l fPdj. This leads to 



• ' ' \(d t + d z )d t h-l&KGT m rfpp\ Xl 

Ao 



V z — = 8ttG / d\EV l ± T lp pP + - \(d t + d z )8 t h - lQnGT pi3 p l p<> 



The terms in square brackets are boundary terms. These may be important at decoupling however due to the finite 
thickness of the last scattering surface we expect these fluctuations to be smeared out on the scales of interest and so 
we neglect these terms. For a string source 2 in the light cone gauge (X° + X 3 = t) this yields 



V 2 ± S = -8ttG M J dcrX.V_L^ 2 '(x - X) 



where the quantities are evaluated at t r = t + z — X 3 (a, t r ). The wavenumber k is related to the multipolc moment, 
I, for I > 60 via k 2 ~ 1(1 + l)/(io — tr) 2 ~ 1(1 + l)Ao> wnere is the conformal time today. 
In Fourier space we have 

- k 2 S k (t r ) = i8nGtik A J daX A (a, t^^^ (2) 
where A — 1,2 runs over the transverse coordinates. 

POWER SPECTRUM 
Power Spectrum on small angular scales I > 500 

In order to calculate the power spectrum we find the unequal time correlator for density perturbations formed at 
different light cone crossing times and then we sum these up between the last scattering surface and today to get the 



1 Note that the formula {dt — d z )f{Z) = E 1 ^ generalises for a general photon path to (dt + P^di) = E 1 ^ 

2 The energy momentum tensor in the conformal gauge is = J dtda{X^X v - X' M X' M )<5( 4 ) (x - X) 
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total power spectrum P(k). Denning the unequal time correlator as 

Pfa.ti.ta) = (2^) 2 <5( 2 )(k 1 + k 2 )( ( 5 kl (<i)4 2 (<2)) , (3) 
we can find it by integrating over the string source terms ^ at the given times, 

P(k,t u t 2 ) = (8^G M ) 2 ^ J dodo' [x\oM)X B {o' M)^^^-^'^ t2 ») , (4) 

where A = (27r) 2 <5(0) is a formal area factor. We obtain the power spectrum through the sum from last scattering t\ ss 
to today <o : 

P(k) = S dt i dt ^P{k, h,t 2 ) 

Now we use the following three assumptions (for the network at a given time as in ref. |10j): (i) The string ensemble 
is a Gaussian process, i.e. we can find all the correlation functions in terms of two point correlators, (ii) We have 
reflection and translation invariance of the transverse coordinates, (iii) We have reflection and translation invariance 
of the worldshcet coordinates. This means that for equal time correlators we can write 

Therefore we can write (for equal time correlators) 

(x A (o,t)X B (o\t)) = °— V (o-o',t) 
(X A {o,t)X' B {o',t)) = —Nh{o-o',t) 

I A \ ft AB 

(x' A (o,t)X' B (o',t)) = —T(o-o',t) 

and hence we have 

((X A (o 7 t)~X A (o' 7 t)) 2 ) = do 1 do 2 T(o 1 ~o 2 ,t)=:T(o-o',t) 

J a 1 Jo 1 

((X A (o,t) - X A (o',t))X A (o',t)^ = J doiMi(oi- o,t) =:Tl(o - o' \t) . 

This is an extremely powerful relation since we can use the fact that at small angles, i.e. small distances for correlators 
we have that T(o,t) oc o 2 ,V(o,t) oc o°,U(o,t) cx o 2 . In particular we write (for small angles) 

T(o,t) w to 2 ,V(o,t) Mv 2 ,Tl(o,t) « %o 2 

where t 2 = | — v 2 and £ = £t is the correlation length of the network which scales in (conformal) time, i.e. £ = constant. 
At time t > ti SSl we have the following relation between the wavenumber k and the multipole I 

that is, the correlation length at time t > t/ ss corresponds to a multipole 500(ti ss /t) since the correlation length at last 
scattering corresponds to a multipole of approximately I « 500. Therefore, for angular scales below 500 we cannot 
integrate back to last scattering and instead to ti ss (500/l) > ti ss . 
For unequal times, we have on small angular scales 

t(o,o'MM) = ((X A (o,t 1 )-X A (o',t 2 )) 2 )^f(o-o') 2 + v 2 (t 1 ~t 2 ) 2 
Hence, on small angular scales, 

P(kMM)^e~ k ^^-t,?l± 
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and so we have that P{k,t\,t2) ~ P(k,ti,ti), such that 

P(fc)« jf ° '"dtiP(Mi,ti) 

where we use the renormalised time t/ti ss . 
Now using 

P(k,t,t) = (8^) 2 ^ J da + J da_ (x A (a,t)X B („\ty k -^-^'^) 
where a + = a + a' and er_ = a — a', we approximate the integral as 

where we note that the range of integration (in a) is limited by the length of string in the network. For large wavenum- 
bers, k > 5000, we may use the small angle approximations as outlined above. However for smaller wavenumbers we 
can approximate this by constraining the range of er_ to (— £/2,£/2), since simulations show that in the range the 
approximations made above are reasonably accurate for k > 500. We can make a better approximation motivated 
by the velocity correlator as found numerically in rcf. [14 3 . In particular we approximate the velocity correlator (see 
Figure to be 4 

V(a,t)*v 2 (l-^f)exp(-|a|/|) 




3 The velocity correlator here is in the light cone gauge. However the behaviour in any other temporal gauge should follow a similar 
behaviour. 

4 This approximation obeys the constraint J daV(a) = which is a reflection of conservation of momentum in the network. 
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Figure 1: Plot of the improved approximation of V(a, t)/(v 2 ). The velocity approaches zero at the correlation length at each 
time t and displays an anti-correlation which is a result of momentum conservation in the network. The plot is in units of the 
correlation length. 

Using this approximation we may safely push the range of integration for V(cr, t) to a G (—00, 00). However, using 
this (better) approximation here makes little difference to the total power spectrum and so we do not use it at this 
junction. Nevertheless, this approximation becomes very important for the bispectrum and trispectrum. 

The second term in the above approximation gives a subdominant contribution to the power spectrum at all angular 
scales and therefore we discuss only the first term below. The power spectrum is therefore given by 

ri 1 rto/Uss d , nit/2 r- ( fp \ 

1 w A £ 2tl Ji t \ 4 J 

since the factor C^/A ~ constant and again we use the renormalised time t/ti ss and the renormalised worldsheet 
coordinate a' — cr/to to simplify notation. For the purposes of finding the total power spectrum we approximate the 
error function by a constant - in particular, by unity since it will quickly reach this value if it is not already equal to 
unity at last scattering. Hence, we obtain the small angle power spectrum (Z > 500): 



Power Spectrum on large angular scales, I < 500 

In ref. [11 , it was shown that if we compute the large angle power spectrum using the contribution of strings at 
the last scattering surface then k 2 P(k) cx k 2 or l 2 Ci cx I 2 on lengthscales above the correlation length (I < 500). 
However, we will show in this section that the late time small angle contribution of strings, i.e. the small angle effect 
of cosmic strings as the network evolves between last scattering and today, is the dominant contribution to the low 
I part of the spectrum. We find this by integrating the unequal time correlator for all times between t > ti ss and to 
and finding where the peak of the spectrum from such times is located. Such contributions are, in effect, super horizon 
effects at last scattering which come inside the horizon at later times. Essentially we are tracking the evolution of the 
peak of the spectrum (located at I w 500 at last scattering) in time. Here we simply note that the contributions on 
superhorizon scales at a time t = tn fall off sufficiently rapidly for k < aH that they can be neglected relative to the 
late time contributions t>tn from subhorizon strings. 

Integrating the string sources over all times since last scattering we obtain the power spectrum relevant on large 
angular scales (I < 500), 

ep(k)«P Cl » (8»CriVs£fyln(>) (7) 

This is only strictly valid, as mentioned above, for I > 60 where the flat sky approximation remains good. However, 
it can be extrapolated to lower I to provide useful estimates as can our subsequent results for the bispectrum and 
trispectrum in these regimes. From this expression and the small angle result, we see that the peak of the power 
spectrum lies near I w 500 with a logarithmic deviation from scale-invariance on large angular scales and a steeper 
fall-off on smaller scales. The large angle result is qualitatively in agreement with early numerical work on this 
problem in ref. |15j (see, more recent results in [5]), though we note that quantitatively the rise towards a peak at 
I m 500 is much steeper for two reasons: (i) The matter-radiation transition breaks scale-invariance because the density 
of strings is substantially higher in the radiation era and is therefore falling through recombination and afterwards 
as the string slowly responds to the smooth increase in the expansion rate, (ii) We do not take into account the 
pre-recombination perturbations seeded in the cosmological fluid by cosmic strings which create additional intrinsic 



7 



CMB anisotropics. Since the power spectrum is not our focus here in this analytic work, we leave these quantitative 
issues for further study elsewhere [P2] . 

BISPECTRUM 

Bispectrum on small angular scales 

The calculation above for the power spectrum ^ can be easily extended to the bispectrum by defining the analogous 
unequal time correlators 

(<5 kl (h)S k2 (t 2 )M*3)) - (2^) 2 <5(k! + k 2 + k 3 )s T (k!, k 2 , k 3 , h, t 2 , t 3 ) , 

which are again found by integrating over the string sources 
B T (k 1 ,-k 2 ,-k 3 ,t 1 ,t 2 ,t 3 ) = 

i (8^) 3 -^^^ J da,da 2 da 3 (x A (a 1% h)X B (a 2 , t 2 )X c (a 3 , i 3 )e lk * J 



where a = 1, 2,3. As for the power spectrum, the bispectrum is found by integration of the unequal time bispectrum 
between last scattering and today U € [i; ss , to] ■ A similar calculation to this was carried out in |10j . except that the 
contribution was considered to be dominated by the strings near the last scattering surface i.e. t « ti ss . Just for clarity 
at the outset we present some of our notation. We are concerned with triangles formed by the three wavevectors 
ki, k 2 , k 3 such that ki + k 2 + k 3 = so we define 

fci = |ki|, k 2 = |k 2 | , fc 3 = |k 3 |, Ki2 = ki-k 2 , k 23 = k 2 ■ k 3 , k 13 = ki • k 3 . 

A particular concern will be the area of the triangle (especially for flattened cases where it vanishes) with 

Area A = ^2(k\k 2 2 + k\k\ + fcfifcf) -k\-k\-k\ = \^jk\k\ - k 12 ... 
K\ = 2AreaA/fci , K 2 — 2AreaA/fc2 , K 3 = 2AreaA / 'k 3 . 

We note also that we will temporarily drop the subscript referring to time, which we will treat as dimcnsionlcss with 
the range of t in [l,toA{»s]- We will continue to use k in place of I, but their relation is given by §5§s in the flat sky 
approximation . 

Writing C ABC = X A (a 1 )X B (a 2 )X c (a 3 ) and D = k a .X a , and again using the assumption of a Gaussian process 
we find that 

(C ABC e lD ) = i ((C ABC D) - (xVi)£>) {x B (°2)d) (x c (a 3 )D)) 
Now writing X\ — X{a\), etc we have 

(C ABC D) = (x A X B )(x^ a .^ a ) + (x A X^)(x B k a .X a 
X B X% ) ( X A k a .X a 



S AB r / 

-j-vfa o 2 ) [{x^k? (xf - x B ) + x°k?(x? x£>) 

S AB 

— V(<7i - <r 2 ) [k° n( CTl - a 3 ) + k%IL{a 2 - <t 3 )] 

+ — V(a 1 - a 3 ) [k? Ufa - a 2 ) + k B U(a 3 - a 2 )] 
S BC 

+ — V(a 2 - a 3 ) [k A n{a 2 - ax) + k A U{a 3 - CTl )] 
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and 

(D 2 ) = (k a .X a k b .X b ) 



{{kt{xt x£) + kf(xf x?))(k*(xf x») + ki(x 3 B x? ))) 

[/tl 3 r((7i - CT 3 ) + H 23 T{(T 2 - 03 ) + Ki 2 r(<7i - <7 2 )] 



2 

where Ki 2 = k 1 .k 2 , etc. We also find that 

fcffc 2 % c (xtD) (X 2 B D) (X°D) = 

^(«12lli2 + Kl3ni 3 )(Kl2ni2 + K 2 3n23)(Ki 3 ni3 + K23II23) 
o 

(where we use the notation II, j — H (<7j — <jj)). The second term is subdominant (as in the case of the power spectrum) 
and so we neglect it in the following. 

In order to calculate this we note that only two of 012, 013, 023 are independent. The choice of independent variables 
depend on the quantities under consideration, e.g. for the term V13II12 we use (712,013 and any term involving 023 
is then written in terms of these variables. Therefore the prescription is to make a transformation from a\, (72,03 
to these variables. Clearly one of the integrations is then independent of the two independent variables chosen and 
integrates to give C, as in the case of the power spectrum. In the following we consider the term V13II12, i.e. 

J da 1 da 2 da 3 V 13 U 12 exp(-D 2 /2) = ^ J da 12 da 13 V 13 Il 12 cxp(-D 2 /2). 

We can assume the small angle approximations for T, II but as we will outline we need to use the better approximation 
in the case of V, as detailed in Section 3. In order to simplify the calculation we note that 

-D 2 t 

— ^— ~ (K12012 + K 130"l3 + K 23(c r l2 - CT13) 2 ) 



t ( ,, K 23_ ^2 , ^2^3 — K 23 2 

4~ 2 12 k^ + ~ 2 °* 



This indicates that we can make the integration separable using 

dn il (]<T V ,Y i3 I\ i ,exp(-D 2 /2) 

t 2 *%°12\ /" j t/ ( 7 2 Arca A 2 



J do 12II12 exp k ^ 12 j J da 13 V 13 cxp | -t 2 

)/■ 



fc 2 



2~ CT 13 



da 12 n 12 exp (^J^tk^ J do 13 V 13 ex P ^-t^Mlj 

where AreaA — \Jk\k\ — k 23 /2 and (for simplicity of notation) we denote Ki = 2AreaA/fci- Since Ki (divided by t ) 
may be less than 500, i.e. correspond to an angular scale below that of the correlation length then we should use a 
better approximation to V\ 2 . From earlier considerations we find that 

J dcrll(a)cxp ^— - = ~^ J da a 2 exp ^— - 

c 2Vierf(fc|t/4) - k£iexp{-{k£i/4:) 2 ) _ c 1 ~. 
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and also 

y daV(a)exp I 1 w 2u 2 y da - exp cr/^J exp I 



„ q-2 V^(2 + i 2 KH 2 )erk(l/(jKt)) exp(l/(^)) 2 - 2£ift _ 2z; 2 

Using this notation (i.e. the above definitions of fi and f 2 ) the total bispectrum on small angular scales (for which 
all ki > 500), after some algebra, reads 

(Kn>4h(k 3 t)M K 3() + Ki3*a/l(*aD/»(tfa£) + la3*?/l(*iD/a(*i{)) (10) 

In order to carry out the time integration we identify ki , Ki as angular multipoles and £ = £t (where we use the 
renormalised time as usual). We have calculated the bispectrum numerically over the full range of multipoles for 
which the flat sky approximation is valid and the results are shown in Figures 3 and 4. The end result is a fairly 
featureless flat bispectrum with very localised and modest upturns for flattened triangles and suppression of squeezed 
triangles because of causality. This is quite unlike the inflationary bispectra found in ref. [5] and should be easily 
distinguishable given a sufficiently significant signal. In the analytic spirit of this paper, however, we press on to give 
some simple analytic approximations to the full bispectrum (??) in all the different regimes: 

Approximation (1): In the small angle limit at high multipoles with kit;, Kit; > 10 at ti ss we have fi(kit) ~ 1\pK 
and h{Kii) » ^l 2 K 2 t . This then implies that 



4 



fto/tlss 

B(k 1 ,k 2) k 3 ) fa - J dt(8irGn) 



- Cj CQV 2 TT k\ + k 2 2 + fcg 

A 64i 4 £ 2 fc 2 fcffc|AreaA 
(S ^ /()! ^ kl + k 2 + kl 



A 64t 4 | 2 fc 2 /cp|AreaA 
3 ££c u 2 7r k\ + kl + kl 1 



The difference between this result and that of ref. f |10j) is due to the approximation made to make the integration 
over the a coordinates separable. This has little effect on the quantitative result and we observe the same asymptotic 
limit. 

Approximation (2): For Kt; > 5 (at the initial time t start) & n d any value of k we find that 

{klklklf^B^MM « (8-G M ) 3 §^J 

bit t 

5( K23 erf(fc 1 ^/4) + Kl3 erf(fc 2 ^/4) + K 12 erf(fc 3 &/4)) 



3(fc 2 fc 2 fc 2 ) 1 /3Area 



A 



where k rn — min(fcj,500) tells us the earliest time, t start , to which we can integrate back and for which k m t 

(k m /500)t/ti ss = 1 at t = t s tart- 

Approximation (3): For the intermediate range Kt; G (3, 5) we find 
{klklklfl 3 B{k x ,k 2 ,k,) * (8^) 3 ^^ 

•A 64t t; 



VbiKwertfatt/Qy/ Kit + K 13 erf(k 2 tt/4,)J K 2 t + K 12 evf(k 3 tt/A)J K 3 t) 



3(fc 2 fc 2 fc 2 ) 1 /3Area 



A 
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This approximation is adjusted slightly for K£, S (1,3) where we replace erf(fc 3 £i/4) by 0.8erf(/c 3 £i/2.5). However 
since we are interested in only a quantitative approximation to the level of signal from the bispectrum we note that 
we can use the expression for Kt; £ (3, 5) in the full range (1,5). 

Approximation (4): For flattened triangles K£ < 1 the behaviour of the bispectrum gives (for small angles fc£ > 4) 
we need an improved estimate finding 



(k\klkl) 2 ^B(k u k 2 ,k s ) » (8ttG M ) 



, c v 2 ir 



A 64t 4 | 2 

V / 5(K 23 erf(fc 1 |i/4) J ft: 1 ^ + K 13 erf(fc 2 £/4)tf 2 | + Kl2 erf(fc 3 ^/4)if 3 f ) 



6(fc2/c2fc2)l/3 AreaA 



(8ttG/i) 



3 ££ c i; 2 7r K 23 erf(fci|i/ 4 ) + Ki 3 erf(fc 2 £i/4) + K 12 erf(fc 3 £i/4) 



A 32i 4 £ 



3(fcffc 2 2 fcl)V3 



This approximation is further motivated by the asymptotic expansion of erfc(x) » exp(— x 2 ) + . . . for x large. For 

fc£ < 4 we note that the formula behaves more like Approx. (3) for lower and lower values of Kt; before it becomes 
accurate. It is important from these approximations to note that although the bispectrum grows as we approach 
collapsed configurations (for which AreaA — > 0) it is not divergent and attains a finite value on the boundary. 

To summarise, let us provide a final simple expression which is valid over all multipoles / < 2000 (e.g. relevant for 
both the Planck and WMAP experiments): 



(fc 2 fc 2 2 fc 3 2 ) 2 / 3 B(A :i ,fc 2 ,fc 3 ) a (SttGm)' 



;££ C Q V 2 1T 



A 64i 4 | 2 

5(K2 3 erf(fci£/4)£>(#i£) + Ki 3 erf(fc 2 £i/4).D(^ 2 f) + K 12 eri(fc 3 £ / f 4)D(K 3 £)) 



3(fc 2 fc 2 fc 2 ) 1 /3Area A 



where the function D is defined by 



D(K£) = 1 



if lf£ > 5 

if e (o, 5) 



(ii) 

(12) 

(13) 
(14) 



Although this expression will not accurately describe the very localised upturn near the edges of tetrahedron, this 
region makes a negligible integrated contribution to Jnl- Going further out beyond I > 2000, these upturned edge 



regions become more significant, so we must replace the second line above (14 1 by the following 



D(KZ) = y/Kt/5 

= k£/VE 



if K£ E (1,5) 
if iff < 1.. 



Bispectrum on large angular scales 

If min ki — k rn < 500, then we must restrict the range of times over which we sum to (t start, to). We find that, 
as in the case for the power spectrum, the bispectrum decreases as a logarithm as we restrict the range, i.e. the 
bispectrum oc ln((fc m /500)to/iz ss ) ~ hi(fc m /10). Therefore, we can include this effect by multiplying the formulae 
in the previous section by ln(fc m /10)/ln(50) for k m < 500. The extension to large angular scales allows us to make 
approximations to the level of non-Gaussianity in the parameter range of Planck. In Figures 5 and 6 we illustrate 
k — fci + fc 2 + fc 3 = const, slices through the bispectrum to reveal the cross-sectional shape using the approximations 
above. 
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Estimator of /jvl 

As detailed in [TB] the signal to noise ratio of the bispectrum is given by 



where &/ 1 ; 2 ; 3 corresponds to the reduced bispectrum, which is precisely the bispectrum as calculated in the flat sky 
limit and where we assume that the noise is cosmic variance dominated. From the images shown in the last section 
we see that in the range of interest for Planck we can find a reasonably accurate estimate for /jv ^ by summing over 
the equal I values. This gives 

2 \ ,a f I I I \ 2 b h,h,h 



(S/Nf * ±£z 



7T 

I 



o o o; c^CfcC, 



1 y, l(l 4 bui) 2 



2yfafi ^ (/ 2 C ; /(2^))3 

For the local model of inflation the signal to noise ratio is proportional to /jvl- Since [5] the local bispectrum is 
proportional to for I < 1000 we normalise the signal to noise ratio of cosmic strings to that of local models of 
inflation, for which 1% U w 2 x 10" 17 7jvl- We have computed / 2 C;/(2tt) using CMBFAST with WMAP5 values for 
the parameters. Therefore, defining, 

j2 _ /^jy^2(cosmic strings) j^gj^ 2 (inflation) 

we plot /jvl against / for multipoles summed between 60 and I. The fluctuations in the plot arise since cosmic strings 
do not strictly follow a local model and because of the transfer functions appearing in the C;s. However, we find that 
for the following values of the cosmic string model parameters, as found in simulations ([T7],[T3], [T5] ) 

v 2 = 0.18; t = 0.42; -f- = 2 

A 

and G/i = 7 x 10 -7 that Jnl ~ — 75cq. Preliminary estimates of Co from numerical simultaions suggest that it is 
significantly below unity, so we estimate Jnl~ —20 (see Figure [2]). Given the most recent estimate of Jnl for local 
models /nl = 38 ± 21 ([19] ). It should also be stressed that cosmic strings produce a very different shape to local 
non-Gaussianity, so analysis of compatibility with WMAP or Planck needs to be specifically investigated ([12J). 



Numerical results 



We have plotted the three dimensional bispectrum in the multipole range of Planck < 2000 by numerically 
integrating the analytic result (see Figures [3] and [4] ) . In particular we plot the quantity (hl2h) 4 ^ 3 B(h,l2,l3)- We do 
not integrate out fully to the edge since the result begins to diverge there and is more difficult to evaluate numerically. 
We find that the signal peaks at around (500, 500, 500) and drops towards the corners. Therefore, our assumption 
that using the equal I values to estimate the level of non-Gaussianity seems valid. However, we do expect that the 
signal rises very near the edge (before it levels out) and so we use our approximations as detailed earlier to track this 
behaviour on slices across the bispectrum. In particular we show the slice for h + h + h — 2000 (Figure [5J and for 
h + h + h = 1500 (Figure [6JI. In order to see the behaviour more clearly we have multplied the bispectrum by a 
negative sign and we have shown the former plot as a contour lot. The slices are plotted in units of e 3 = (87rG/i) 3 . 
There is a peak towards the edge but only a factor of a few that at the centre point. Also it begins to diverge very 
near the edge for these slices. Therefore, we feel justified in neglecting the edges in our estimator for /jvl- 
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Plot of f , against I 
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Figure 2: Measure of the level of non-Gaussianity using f^ al 

TRISPECTRUM 
Trispectrum on small angular scales 

As in the section on the bispectrum we regard k as the angular multipole and £ to be the correlation length as 
measured at the original time (normalised by some factor as discussed previously). We will return the to the time 
dependence later in the calculation. However, we note here that as in the case for the power spectrum and bispectrum 
the total trispectrum reduces to the sum of equal time four point correlators between t start > ti ss and to- 

The trispectrum, in the flat sky approximation, is given (at any given time) by 

AAA 3 <*kJ c - (2 7 r) 2 ^ 2 )(k 1 +k 2 +k 3 + k 4 )T c (k 1 ,k 2 ,k 3 ,k 4 ) 

where the subscript c denotes the connected part. This suggest that we normalise the trispectrum by the area factor A. 
However, as we shall see below we should in fact normalise by A 2 for parallelogram configurations. The unconnected 
part, denoted hereafter by the subscript uc, is given by 

AAA 3 <Sk 4 ) MC = (2tt) V 2 )(k! + k 2 )^(k 3 + k 4 )P(fc 1 )P(fc 3 ) + (2vr) V 2 )(k x + k 3 )5^(k 2 + k 4 )P(fc 1 )P(fc 2 ) 
+(2n) 4 S^(k 1 + k 4 )<5^(k 2 + k 3 )P(fc!)P(fc 2 ) 
= ^l 2 A.fcAa.fc^i)^) + S kuk3 S k2M P(k 1 )P(k 2 ) + S kuk3 S k2M P(k 1 )P{k 2 )) 

We will calculate the four point correlator and subtract off this unconnected part to find the trispectrum. 

In this paper we will consider only parallelogram configurations. In this regard we are motivated by the earlier 
section on the bispectrum which showed that as the triangle collapses it grows (though does not become infinite). 





Figure 4: Plot of the 3D bispectrum (hhh) 4 ^ 3 B(h, h, h) with difference slicings to that of Figure^ 



14 




Figure 5: Contour plot of the slice through the bispectrum at £\ U — 2000. 



Bispectrum for 1 +1 +1 =1500 
12 3 




Figure 6: Plot of the slice through the bispectrum at U — 1500. 

Clearly, the parallelogram shape should replicate this behaviour. 

For planar parallelograms (planar shapes are, of course, appropriate since we are using the flat sky approximation) 
we have that opposite sides are antiparallel. To impose this we need to use two delta functions and this implies that 
the trispectrum is normalised by A 2 . This is true in general, however, which becomes more apparent also when we 
note that we can decompose the delta via 

^(k! + k 2 + k 3 + k 4 ) = / d 2 L<S (2) (ki + k 2 + L)$W(ki + k 2 - L). 



Then we infer that we can approximate the total trispectrum by summing over the possible parallelogram configura- 
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tions (to specify a parallelogram configuration it suffices to have two of the sides and a diagonal) 5 . 
Motivated by the above discussion we define 

r T (k 1 ,k 2 ,k 3 ,k 4 ) = (SkJkJkJkJ 

where the multipoles trace out a parallelogram configuration. 

In a similar fashion to the calculation of the power spectrum and bispectrum, we find that 



T T (k l7 k 2 ,k3,k 4 ) = 



Writing C ABCD = X A (<Ti)X B '(a 2 )X c (a 3 )X D '(0-4) and E = k a .X a and using again the assumption of a Gaussian 
process we get that 

l^jABCD e iE^ _ l^jABCD _ \qABCD j^2 _|_ }_(jABCD g4 _|_ ^ 



Now 



, . S.AB zCD zAC XBD zAD sBC 

(c ABCD ) = (x A x B x§x?) = d —v 12 d —v 34 + d —vJ—v 24 + 6 —vJ—v, 



2.3 



=► k A k B k^k B (C ABCD ) = ^(k 12 k 34 V 12 V 34 + k 13 k 24 V 13 V 24 + k 14 k 23 V 14 V 23 ) 

Next we consider (writing C ABCD as ABCD) 

(ABODE 2 ) = (C ABCD ) (E 2 ) + 2 (AB) (CE) (DE) + 

where denotes the contribution from terms found by permutating the symbols A, B, C, D in the the second term 
on the right hand side of the equality 6 . Now we note that (using J2 i k^ = 0) 

(AB) (CE) (DE) = <^ tx° £ kf Xf\ tx° £ kfxf\ 



2 

5 AB Vi 



J¥3 / \ 
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where we use the fact that 11^ = (xft-X^j = — < X A X A t >. Suppressing the summation notation we find that 
k A k B k^kf (ABODE 2 ) = (C ABCD ) (E 2 ) + 

^ (Vl 2 K 12 K 3j Il 3j K 4j Il 4j + V 13 K 13 K2jU2jK 4j U 4j + Vl 4 K 14 K 2j Il 2j K 3j Il 3j ) + 
^ (V23 K 23 K lj^lj K 4j^4j + V24K24KljN-ijKsjIlsj + V34/t34/tijIIijK;3 :? n3j) . 

Now we proceed to investigate 

(ABODE 4 ) = (ABODE 2 ) (E 2 ) + 24 (AE) (BE) (CE) (DE) . 



We shall see later in the section that the Gaussian contribution is given by square configurations. 
We note that the factor of two appears since (AB) = (BA). 



1G 



Wc have seen above that (AE) = — ^ J2j^i kfH-ij so (again suppressing the summation notation) we have 

24 

fc^fcf fcf fcf (AE) (C£) (Df?) = — K lj U lj K 2j Il 2j K 3j U 3j K 4j U 4j 
Finally, we collect the terms to find that 

fc^fcf fcf fcf (C ABCD e iB ) - [1 (K12K34V12V34 + K13K24^1 3 V 2 4K14K23^14V 23 ) 



where we have 

(e 2 ) = ({ktxA + kixA+kfxAxkZxB + kixB + kfx*)) 
= ~2 («i2ri2 + K i3 ri 3 + Ki 4 ri4 + K 23 r 23 + K 24 r 24 + K 3 4r 34 ) 

and where we use the notation to represent the non-equivalent contributions by permuting the symbols (1,2,3,4) 
of the term in the brackets. The dominant term on all angular scales is found to be given by the first term (~ 
qABCD e -E / 2 -j anc ^ therefore, in what follows we neglect the remaining terms. 

For a parallelogram configuration we have K\ 3 = — fc 2 , k 24 = — fc 2 , k 23 = k 14 = — k 34 = — k 12 . This implies that, in 
the small angle approximation, 

(E 2 ) = — (ki 2 CT 2 2 - fc 2 r i3 - Ki 2 (T 2 4 - Kl 2 CT 23 - fc 2 CT 2 4 + K 12^3i) 

Again the choice of independent variables is dependent on the integrand under consideration. To elucidate this we 
consider the term 

J daida2da 3 da i Vi 3 V24 : exp(-E' 2 /2) = — J c£ai 2 cfoi 3 <&7 2 4Vi 3 IIi 2 exp(-E 2 /2) 

In a similar manner to the case of the bispectrum we can simplify the calculation by making a separable approximation. 
In particular we find 

exp(- (E 2 ) /2) = exp ^— (ki 2 (t 2 2 - k 2 Ti 3 - Ki 2 cr 2 4 - ki 2 ct 23 - k\a\ 4 + Ki 2 cr 2 4 )^ 



f t ( . Kyi , 2 . ^1^2 K 12 2 \ 
ex P ( \ (^1^13 - -JT^) + ^2 U 2i) 



f C 2 f ( k 2 t z cr 2 \ f ' ' " 2 + 2 - 2 

/ daida2da 3 da4,Vi 3 V24, exp(— E 2 /2) = — doi 3 Vi 3 exp I 1 - 13 J / dcr^V^exp 



K{t a 2 4 



where Ki = 2Area,/\/ki as defined in the section on the bispectrum (note that the area of the triangle is half that of 
the parallelogram). Since K 2 (divided by to) may be less than 500, i.e. correspond to an angular scale below that 
of the correlation length then we should use a better approximation to V 2 4. Following these approximations for all 
terms and using the results of the previous sections we find that 

T T (k 1; k 2 ,k 3 ,k 4 ) « |2fc^|fc2 Jq (/ ^i2^ 34 2«: 2 2 Fi 2 y 3 4Cxp(- J B 2 /2) + J da 13 da 24 k 2 k 2 2 V 13 V 24 exp(-E 2 /2)^j 

We note also that in the case of a square this gives P(ki)P(k 2 ) at each time. This is the Gaussian contribution. 
Therefore the approximation to the total trispectrum is given by the sum over non-square parallelogram configurations. 
In order to compute the integral we use the more accurate approximation for V. Using the previous expressions for 
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the separate integrals we find 

rr( kl ,k 2 ,k 3 ,k 4 )^ f ^£j l _ 1 (242h{K 2 hh{k 2 i) + k\klf 2 {K l 0f 2 {k l i) 

A 2 gk\k\k\k\ 4 f| 4 Area 3 v 

where Area = 2AreaA- 

This result must be integrated as usual to give the total trispectrum due to cosmic strings between last scattering 
(or if k rn < 500 between t star t = ^ss500/fc m ) and today. Returning the time dependence using the renormalised 
time t/ t s tart and understanding the multipoles to refer to angular multipoles we have that the total trispectrum 
contribution by the parallelogram considered here is 



rto/t 



C 2 e (8ttG M ) 4 v 4 1 
A 2 WkJkJkJkJ^ffAxea 3 



In order to understand the behaviour for the trispectrum we note the following approximation to the integral given 

by 

1 ' A 2 ek\klk\k\ 4 J 2 Area 

(2/t? erfc ( t= — | exp ( — | erfc ( ^ | cxp ( — 

+fc 2 fc 2 erfc f -J— ) cxp [ — | erfc | — | exp [ — — ^ )). 

" UkiiJ V(£M) 2 / VM/ V(^i<) 2 / 

Using the asymptotic expansion of erfc we find that as K — > the formula is independent of K, i.e. of Area. Therefore, 
as for the bispectrum, the trispectrum grows towards the edge but is not divergent (i.e. it reaches a cutoff). 



Trispectrum on large angular scales 

If min ki = k m < 500, then we must restrict the range of times over which we sum to {t sta rti to)- Again we find that 
the trispectrum decreases as a logarithm as we restrict the range, i.e. the trispectrum cx ln((fc m /5QO)to/i; S s) ~ 
ln(fc m /10). Therefore, we can include this effect by multiplying the formulae in the previous section by 
ln(fc m /10)/ln(50) for k m < 500. This is the same effect as was evident in the case of the bispectrum and the 
power spectrum. That we can extend the range of multipoles to less than 500 means that we can use this formalism 
to make a prediction for tjvx generated by cosmic strings in the multipole range of Planck. 



Estimator of tnl 

The signal to noise ratio for the trispectrum (assuming the noise is cosmic variance dominated) is given by 

\T h i 2h i 4 (L)\ 2 1 ^ \ T hi2hh( L )\ 2 



(S/NY 



h <ih<UL (^ + l)C h C l2 C h C h 24^2. { 2L + l)C h C l2 C h C h 



where T; 1 i 2 / 3 ; 4 (L) is one of the possible configurations of the full-sky trispectrum. In order to find the signal to noise 
ratio we therefore must relate the full sky trispectrum to the flat sky equivlaent (|20j). 
Decomposing the temperature perturbation in terms of spherical harmonics we have 

AT 

— (n) = ai m Y lm (h). 
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In the flat sky limit the temperature perturbation is expressed in terms of Fourier harmonics 

AT f d 2 l , a& 

The relationship between the full-sky and flat-sky limit is obtained by noting 



5(1 1 + I 2 + 1 3 + I 4 ) = J ^_ e ^+i 2 +i3 + U). fl 

m 

(a hmi a l2m2 a hm3 a hmi ) c = (jlf =1 ^ e - im ^ (2K) 2 5(h + h+l 3 + h)T^$ 

where </>; denotes the polar angle of 1 and T„*'f*} denotes the flat sky trispectrum given by a configuration with 
multipoles h,l2,ht '4 (and where 1 12 is the diagonal of the quadrilateral that makes a triangle with Using these 

formulae we find 

<K(li + b + 1 3 + U)) = / d 2 LS((h + h + L))*((Ia + 1 4 - L)) 

( 2 ^?S( ntl v w J ** 

V(2/i + l)(2/ 2 + l)(2/ 3 + l)(2/ 4 + 1) 

(—i) M 1 ^ ^ 2 ^ 1 ( ^ 3 ^ 4 ^ 

I mi fni M I \ TO3 TO4 — M 
Combining these results we can write the full sky trispectrum in terms of the flat sky trispectrum as 



h h L \( h k L 




T, 



W = ^ 1 V(^ + l)(2/ 2 + l)(% + l)(2/ 4 + l)( J} ^ J ) ( '» £ L Q ) T£g(L). (15) 



With this consideration we can make an estimate of the signal to noise ratio for the trispectrum. We note also that 
we need only consider even terms since the Wigner-3j vanishes for ij + lj + = odd. 

We are interested in finding an estimate for t^l- We note that for the local model of inflation, [ST] (for which 
tnl = (Q/nl/S) 2 ), the signal to noise ratio is given approximately by 

{S/N) local M 5 x ^^looal)^ ^ ^ % ^S^local) ^ 

Observationally the trispectrum has not received much attention with the rather weak constraint tjvl < 10 8 [? ]. 
Recently there has been a significant improvement claimed using N-point probability distributions |22j with — 5.6 x 
10 5 < tnl < 6.4 x 10 5 . Again these results are for a local-type non-Gaussianity which is unlike the less peaked cosmic 
string trispectrum, so a more specific analysis will be necessary to achieve quantitative constraints. 



CONCLUSIONS 



In this paper, we have endeavoured to analytically calculate the cosmic string power spectrum, bispectrum and 
trispectrum over a range of scales relevant for CMB experiments on both large and small angular scales. We have been 
particularly focused on extending previous work to multipole ranges (/ < 2000) applicable for the Planck satellite, an 
experiment which has the potential to dramatically improve constraints on all these correlators. We have presented a 
relatively featureless shape for the bispectrum over the relevant range which should be easily distinguishable from the 
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oscillatory peaks of inflationary bispectra [H] if there is a significant signal discovered. Our preliminary estimates of 
/atl from cosmic strings indicate that Planck constraints from the bispectrum should be competitive with those from 
the power spectrum. We note that we obtain a considerably smaller estimate than the Jml — —1000 in ref. [11] for 
several reasons, including the extension of our analysis over lower multipolc ranges relevant for WMAP, a more careful 
comparison with /jyi, estimators used in the literature, and a lower normalisation of the string spectrum derived for 
Nambu strings [5], rather than that obtained from field theory simulations [BJ. Nevertheless, as we shall discuss, our 
estimate contains many uncertainties and more detailed and accurate forecasts are the subject of ongoing work |12j . 

We have also evaluated the CMB trispectrum (for parallelogram configurations) on the relevant multipole scales 
for WMAP and Planck. Again we find a relatively constant trispectrum, finite in all regimes and for which we 
do not expect dramatic features for other configurations. Our results indicate a relatively larger signature for the 
trispectrum, in contrast to the amplitude of the bispectrum which is suppressed because it depends on the poor 
correlation between the string velocity and curvature. Our preliminary estimate of tml a few times 10 4 needs a more 
detailed analysis to characterise key uncertainties more carefully. Nevertheless, the trispectrum deserves much closer 
scrutiny observationally and the prospect of constraining cosmic strings should motivate the development of suitable 
estimators. 

These analytic calculations of the post-recombination gravitational effects of cosmic strings offer important physical 
insights into the CMB correlations they induce. However, we note that there are many directions in which they can 
be substantially improved. Detailed numerical investigations of the CMB power spectrum created by cosmic strings 
and other topological defects already takes into account a far wider range of physical effects, notably recombination 
physics around decoupling and a better description of the evolving string network. It should be possible to similarly 
develop the unequal time correlation methods presented for the late time GKS signatures here to calculate both the 
bispectrum and the trispectrum to high accuracy numerically. In the meantime, however, it is a matter of comparing 
the present analytic results with CMB maps induced by cosmic string networks on both large (full sky) and small 
scales in order to get a more accurate normalisation and characterisation of the bispectrum and trispectrum. Given 
the stark contrast between the string shapes with those predicted by inflation, there needs also to be a specific search 
for these signatures in present and forthcoming CMB data, a project which is being actively investigated [T^] . There 
seem to be good prospects of using forthcoming CMB data to obtain new insight into cosmic string scenarios using 
higher order correlators. 
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